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ABSTRACT 



We calculate the nuclear modification factor (Raa) of light hadrons by taking into account 
the initial state momentum anisotropy of the quark gluon plasma (QGP) expected to be 
formed in relativistic heavy ion collisions. Such an anisotropy can result from the initial 
rapid longitudinal expansion of the matter. A phenomenological model for the space time 
evolution of the anisotropic QGP is used to obtain the time dependence of the anisotropy 
parameter ^ and the hard momentum scale, phard- The result is then compared with the 
PHENIX experimental data to constrain the isotropization time scale, riso for fixed initial 
conditions (FIG). It is shown that the extracted value of Tiso lies in the range 0.5 < riso < 1.5. 
However, using fixed final multiplicity (FFM) condition does not lead to any firm conclusion 
about the extraction of the isotropization time. The present calculation is also extended to 
contrast with the recent measurement of nuclear modification factor by ALIGE collaboration 
at ^/s = 2.76 TeV. It is argued that in the present approach, the extraction of Ti^o at this 
energy is uncertain and, therfore, refinement of the model is necessary. The sensitivity of the 
results on the initial conditions has been discussed. We also present the nuclear modification 
factor at LHG energies with ^/s = 5.5 TeV. 



1 Introduction 

The partonic energy loss in a QGD plasma has received significant attention in recent years. 
Experimentally, the partonic energy loss can be probed by measuring the high pt hadrons 
emanating from ultra-relativistic heavy ion collisions. This idea was first proposed by Bjorken 
[1] where 'ionization loss' of the quarks and gluons in a QGD plasma was estimated. In fact, 
the 'stopping power' (dE/dx) of the plasma is proportional to y/e, where, e is the energy 
density of the partonic medium. Therefore, by measuring various high px observables one 
can probe the initial parton density [1]. 

Hard partons, injected into hot QGD medium, can dissipate energy in two ways, viz., 
by two body collisions or via the bremsstrahlung emission of gluons, commonly referred to 
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as coUisional and radiative loss respectively. For electromagnetic processes, it is well known 
that at large energies, radiative losses are much higher than the collisional loss. 

The possibility of QGP formation at RHIC experiment, with initial density of 5 GeV / fm'^ 
is supported by the observation of high pt hadron suppression (jet-quenching) in the central 
Au-Au coUisions compared to the binary-scaled hadron-hadron collisions [2] . The phenomena 
of jet-quenching has been investigated by various authors [2]. Apart form jet-quenching, 
several possible probes have been studied in order to characterize the properties of QGP. 

However, many properties of QGP are still poorly understood. The most debated ques- 
tion is whether the matter formed in the relativistic heavy ion collisions is in thermal equi- 
librium or not. The measurement of elliptic flow parameter and its theoretical explanation 
suggest that the matter quickly comes into thermal equilibrium (with Ttherm < 1 fm/c, where 
Ttherm. is the time of thermalization) [3]. As for example, one of the major difficulty is to mea- 
sure the thermalization {Ttherm) and isotropization (rj^o) time of the QGP. On the one hand, 
the success of ideal hydrodynamic fits to experimental data [3] implies rapid thermalization 
of the bulk matter created at RHIC. On the contrary, perturbative estimation suggests rel- 
atively slower thermahzation of QGP [4]. However, recent hydrodynamical studies [5] have 
shown that due to the poor knowledge of the initial conditions there is a sizable amount 
of uncertainty in the estimate of thermalization or isotropization time. It is suggested that 
(momentum) anisotropy driven plasma instabilities may speed up the process of isotropiza- 
tion [6], in that case one is allowed to use hydrodynamics for the evolution of the matter. 
However, instability-driven isotropization is not yet proved at RHIC and LHC energies. 

In absence of a theoretical proof favoring the rapid thermalization and the uncertainties 
in the hydrodynamical fits of experimental data, it is very hard to assume hydrodynamical 
behavior of the system from the very beginning. Therefore, it has been suggested to look 
for some observables which are sensitive to the early time after the coUision. For example, 
jet-quenching vis-a-vis energy loss of partons could be an observable where the initial state 
momentum anisotropy can play important role. It is to be noted that the calculations of 
energy loss in Ref. [2] have been performed in isotropic QGP which is true immediately after 
its formation [7]. However, subsequent rapid expansion of the matter along the beam direc- 
tion causes faster cooling in the longitudinal direction than in the transverse direction [4]. 
As a result, the system becomes anisotropic with (pl^) << (pt^) in the local rest frame. At 
some later time when the effect of parton interaction rate overcomes the plasma expansion 
rate, the system returns to the isotropic state again and remains isotropic for the rest of the 
period. Thus, during the early stage the plasma remains anisotropic and any calculation of 
energy loss should, in principle, include this aspect. The collisional energy loss in anisotropic 
media for heavy fermions has been calculated in Refs. [8, 9]. In these calculations it is found 
that the deviations from the isotropic results are of the order of 10% for ^ = 1 (^ is the 
anisotropy parameter) and of the order of 20% for C, = 10. It is observed that the collisional 
energy loss varies with the angle of propagation by upto 50%. 

Radiative energy loss in anisotropic QGP has recently been calculated in Ref. [10] in first 
order opacity expansion where the scatterers are treated as static. It should be noted that 
the energy loss of a parton in anisotropic media depends both on the anisotropy parameter 
and the direction of propagation with respect to the anisotropy axis n. In that case, general 
expression for the two-body potential should be used [10]. Few comments about the calcu- 
lation of radiative energy loss are in order here. First of all we have considered scattering 
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from static charges in which case only longitudinal gauge bosons are exchanged rendering 
the screened potential stable. The effect of anisotropy would be significantly different if one 
considers the recoil of the scattering centre. In case of moving scattcrcr, transverse (mag- 
netic) gluon exchange results plasma instabilities which has to be taken into account. It is 
worthwhile to mention that the growth of unstable modes in anisotropic media greatly in- 
fluences the transport co-efficients and hence the energy loss. This is recently demonstrated 
in Rcf. [11]. Considering a two-stream plasma the authors in [11] show that the momentum 
broadening grows exponentially with time as the spontaneously growing fields exert an ex- 
ponentially growing influence on the propagating parton. In an evolving plasma this aspect 
is an important component without which the results remain somewhat like a zeroth order 
approximation. However, in the present work, we do not include this effect and assume that 
the fragmenting partons propagate in the direction of anisotropy. 

For a parton propagating in the direction of anisotropy it is found that that the fractional 
energy loss increases by a factor of 1.5 - 2 depending upon the anisotropy parameter ^. In 
the present work, we shall apply this formalism to calculate the nuclear modiflcation factor 
for light hadrons. A phenomenological model for the space-time evolution will be used for 
the time evolution of ^ and phard- Since the role of coUisional energy loss in the context of 
RHIC data is not settled yet, we shall not include this process in the present work. 

The plan of the paper is the following. In section 2 we briefly recall the necessary 
ingredients to calculate radiative energy loss in anisotropic media. Then we discuss how this 
can be implemented to calculate ii!AA along with space-time model for anisotropic media. 
Section 3 will be devoted to discuss the results. Finally, we conclude in section 4. 

2 Formalism 

2.1 Radiative energy loss 

In this section, we briefly mention the formalism of the radiative energy loss in an inflnitely 
extended anisotropic plasma (see Rcf. [10] for further details). 

We assume that an on-shell quark produced in the remote past is propagating through 
an inflnite QCD medium that consists of randomly distributed static scattering centers 
which provides a color-screened Yukawa potential originally developed for the isotropic QCD 
medium given by [12] 



with v{qn) = 4:7ias/{(^ + /x^), where n is the Debye mass. a;„ is the location of the nth 
scattering centre, T (summed over a„) denotes the colour matrices of the parton and the 
scattering centre. It is to be noted that the potential has been derived by using hard thermal 
loop (HTL) propagator in QGP medium. In a plasma with momentum anisotropy the two 
body interaction, as expected, becomes direction dependent. In addition, it depends on the 
anisotropy parameter ^ and the hard momentum scale Phard in the following way: 



Vn = V{qn)e 



27r(5(g°)T;(g„)e-^«"-^"r„„(i?) ^ T„„(n). 



,0 



(1) 



m 



1 



(2) 



gP/Phard-\/l+^(P'")^ ± 1 
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n-fP) - ^ I -^v^ ^ (a^^ - (3) 



where, the parameter ^ is the degree of anisotropy parameter (—1 < ^ < oo) and is given by 

C = (Pt)/(2(P2)) — 1- It is to be noted that ^ can also be related to the shear viscosity [13]. 

To calculate the two body potential appearing in Eq.(l) one starts with the retarded 
gluon self energy expressed as [14] 

d'k dfjk) / , _ v'^P^ 
(27r)3 dKf^ Y P-v + ie^ 

We have adopted the following notation for four vectors: = (po^p) = {po,P,Pz), i- e. p 
(with an explicit vector superscript) describes a three-vector while p denotes the two- vector 
transverse to the z-direction. 

To include the local anisotropy in the plasma, one has to calculate the gluon polarization 
tensor incorporating anisotropic distribution functions of the medium. This subsequently 
can be used to construct HTL corrected gluon propagator which, in general, assumes very 
complicated from. Such an HTL propagator was first derived in [15] in time-axial gauge. 
Similar propagator has also been constructed in [16] to derive the heavy-quark potential in an 
anisotropic plasma, which, as we know, is given by the Fourier transform of the propagator 
in the static limit. 

The self-energy, apart from momentum P^, also depends on a fixed anisotropy vector 
n^{= {l,n)) and 11'^'^ can be cast in a suitable tensorial basis appropriate for anisotropic 
plasma in a co- variant gauge in the following way [16]: 

U^"" = a A^"" + /3 B^'' + -f C^'' + S D^"" (4) 

where the basis tensors are constructed out of p^, and the 4- velocity of the heat bath 
u'^. The detailed expressions for the quantities those appear in Eq.(4) can be found in 
Ref. [15, 16]. The anisotropy enters through the distribution function given earlier (see 
Eq.(2)). 

Since the self-energy is symmetric and transverse, all the components are not independent. 
After change of variables {p' = p^[l + ^{p ■ n)"^]) the spatial components can be written as 

= / ^ ( 6^^ + (5) 

^ J Att l + i{v-nY \ P-v + ie) ^ ^ 

Now a,/3,7 and 6 are determined by the following contractions: 

p'W^p> = 

A'^n^U'^p> = {f-{n-Pf)5 



p2 

Trtf^' = 2q; + ^ + 7 (6) 

where the expressions for a, /3, 7 and S are given in Ref. [15]. 

After knowing the gluon HTL self-energy in anisotropic media the propagator can be 
calculated after some cumbersome algebra [16, 17]: 

4 2 \ 

+Ag[{P' -a- 7)^3'^'' + {co' - ^)C^^ + S'^^D^^] - —.P^P" (7) 
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where 



A-l = (p2 _ ^ _ ^)(^2 ^2[p2 _ . p^i 



(8) 



It is well known that the momentum space potential can be obtained from the static gluon 
propagator in the following way: 



= 9 



(g^ + + ) (g^ + m^) - mj 



(9) 



where, in general, the expressions for m^, m^, and are lengthy [16]. For Qz — (which 
is the case here) these simphfy to [10] 



m 



m. 








1 



(10) 



for a parton propagating in the anisotropy direction. Here /j, is the Debye mass given by 
^2 ^ (7^Phard (1 + ^f/Q): whcrc A^f IS the number of flavours. In such case the two-body 
potential in anisotropic media simplifies to [10] 



with 



^(q,0 



1 



1 



q2 + i?(0/x2 
tan"^ 



+ 



(11) 



(12) 



.1+e ■ 

Now the parton scatters with one of the colour centre with the momentum Q = (0, q, g^) 
and radiates a gluon with momentum K = (w, k, fc^). The method for calculating the 
amplitudes is discussed in Refs. [18, 19, 20] and we shall ,for the sake of brevity, quote the 
main results only. The quark energy loss is calculated by folding the rate of gluon radiation 
{r{E)) with the gluon energy by assuming a; + go ~ ^- In this approximation one finds [18], 



dE E f dV 

— — — — / xdx— 
dL Dji J dx 



(13) 



Here Dr is defined as [ta,t^[tc,ta] = C2{G)CrDr, where C2{G) ^ 3, Dr ^ 3 and [ta,tc] is 
a color commutator (see [18] for details), x is the longitudinal momentum fraction of the 
quark carried away by the emitted gluon. 
Now in anisotropic media we have [10], 



X 



dr CrUs L r d'^k rf^q 

dx TT X J TT TT 



q 



(k + q)2 + x2 k2 + x 



(14) 
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In the last expression, v{q,^) is the two body quark-quark potential given by Eq.(ll) 
and X = i^l^^ + 'ml, where = A*^/2 and = /i^/6. 

For arbitrary ^ the radiative energy loss can be written as [10] 

E TT^ XJ '''''' '^'(q2 + i?(0//2)2 [ 2 

_|_ q^ ~ + X 



kl + X 2^q4 + 2q2(x-A;^) + (A;^ + x)^ 
:ln 



^ (q^ -kl + 3x) + vTTfVq4 + 2q2(x-fc^) + {kl + x) ^ 

(15) 

In the above expression, A denotes the average mean free path of the quark given by 

which in this case would be ^ dependent. Note that A^, and A^ correspond to the contributions 
to the mean free path of the propagating quark coming from q-g and q-q scatterings. 
Explicitly with Eq.(ll) we have, 

. -1 _ CRC2{i)Pi f o Aal 

where Cr — 4/3, C2{i) is the cashimir for dj-dimensional representation and C2{i) = {N^ — 
l)/{2Nc) for quark and C2(i) = Nc for gluon scatterers. dA — — 1 is the dimensionality 
of the adjoint representation and pi is the density of the scatterers. Using — pf°/^/T+^ 
we obtain 

1 ^ 18a,pha.dC(3) 1 1 + Nf/6 , . 

X TT^VT+I i?(0l + ^F/4 ^ ^ 

For ^ — >■ Eq.(18) reduces to well-known results [18] 

1 _ 18a,TC(3) 1 + Nf/6 
X ~ 7r2 1 + Nf/4. 

2.2 Hadronic pr spectrum 

The high pt inclusive hadron spectrum in a heavy ion collision can be calculated in a pQCD- 
improved parton model. There are various approaches about how to incorporate the energy 
loss in the hadron production from jet fragmentation. The differential cross-section for 
hadron production is [21], 

E^{AB^iet + X) = K ^ j dxadxbGa/hA{xa,Q^)Gb/hB{xb,Q^) 



(19) 



abed 



X -^{ab^ cd)S{s + i + u), (20) 
TT dt 
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Expressing the argument of the 5-function in terms of Xa and Xb and doing the Xb integration 
we arrive at the final expression: 

E^{AB ^iet + X) = KY^f dXaGa/hA^a,Q') Gb/hA^b,Q') 



abed 



X - ^"^^ ^(ab^cd) (21) 
TT 2xa - xrey dt 



where Xh = {xaXre'^) /{2xa - xtc^), xt = 2^Tl\fs and Xamin = (xTef)/(2 - xtc"^). It 
should be noted that to obtain single particle inclusive invariant cross-section in hadron- 
hadron collisions, the fragmentation function D^jJyZ^ Q^) must be included. Multiplying the 
result by the nuclear overlap function for a given centrality one can obtain the pt distribution 
of hadrons in A — A collisions. However, the inclusion of jet-quenching as a final state effect in 
nucleus- nucleus collisions , can be implemented in two ways: (i) modifying the fragmentation 
function [22] and (ii) modifying the partonic pt spectra [23] but keeping the fragmentation 
function unchanged. We shall concentrate on the former approach, as it is easier to include 
the final state effect in presence of anisotropy. Now the energy loss of high energy quarks 
and gluons traveling through dense colored plasma can measure the integrated density of 
the colored particles. This uon-Abclian energy loss is a function of parton opacity L/X. We 
use the expression given by Eq.(15) which is derived to first order in opacity. 

Now in order to obtain the hadronic Pt spectrum in A — A collisions we modify the 
fragmentation function to obtain an effective fragmentation function as follows: 



Dh/,{z,Q^)^-Dnf,{z\Q^) (22) 

where, z* = z/{l — AE/E) is the modified momentum fraction. Also in order to take into 
account the jet production geometry we assume that all the jets are not produced at the 
same point and the path length traversed by these partons before fragmentation are not the 
same. It is also assumed that the jet initially produced at (r, </>) leaves the plasma after a 
proper time {ti) or equivalently after traversing a distance L (for light quarks ti ^ L) given 



by L{r, 0) = y — f^'^ sin 0^ — Rt cos (p, where R^ is the transverse dimension of the system. 
Now the hadron pt spectra depends on the path length the initial parton must travel and 
the temperature profile along that path. As mentioned this is not the same for the all jets 
as it depends on the location where the jet is produced. Therefore we have to convolute 
the resulting expression over all transverse positions. Since the number of jets produced at 
r is proportional to the number of binary collisions, the probability is proportional to the 
product of thickness functions : 

V{f^^Ty,{r)TBif^ (23) 

For a hard sphere V{r) is given by 

where / d^rV{r) = 1. Using all these and noting that the path length is not a measurable 
quantity we obtain the pr spectra of hadrons as follows: 
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E / <frV{r) [ 



*i dt 



th — U 



E 



dz 



dN 
d^p-f ' 



(25) 



The quantity is the initial momentum distribution of jets and can be computed using 

LO-pQCD (see Eq.(21)). For the fragmentation function Eq.(22) has been used. We use the 
average value of distance traversed by the partons, (L) given by 



/o^- rdr J,'-Li(l),r)TAAir,b = 0)d(j) 



r27r 



(26) 



where (L) ~ 5.8(6.2)/m for RHIC (LHC). The nuclear modification factor, Raa is defined 



as 



dN 



RAAiPx) 



cPpTdy 



dN 



7rO(i7) 



d^PTdy 



(27) 



where the suffix '0' in the denominator indicates that energy loss has not been considered 
while evaluating the expression. 




Figure 1: (Color online) Fractional energy loss for light quark when mean free path is 
independent of ^ (left panel). The ratio of the fractional energy loss of anisotropic media to 
that in isotropic media is also presented (right panel) 
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2.3 Space-time evolution 

For an expanding plasma the anisotropy parameter ^ and the hard momentum scale Phard 
(appearing in Eq.(15) via A as well as in the expression of the Debye mass) are time depen- 
dent. Thus to calculate Raa one needs to know the time dependence of Phard and We 
shall follow the work of Ref. [24] to evaluate the px distribution of hadrons from the first 
few Fermi of the plasma evolution. Three scenarios of the space-time evolution (as described 
in Ref. [24]) are the following: (i) Tiso = t,, the system evolves hydrodynamically so that 
^ = and Phard can be identified with the temperature (T) of the system (till date all the 
calculations have been performed in this scenario), (ii) Tiso oo, the system never comes to 
equilibrium, (iii) Tiso > and Tiso is finite, one should devise a time evolution model for ^ and 
Phard which smoothly interpolates between pre-equilibrium anisotropy and hydrodynamics. 
We shall follow scenario (iii) (see Ref. [24] for details) in which case the time dependence of 
the anisotropy parameter ^ is given by 



e(r,<5) = (- 



(28) 



where the exponent 6 = 2 (2/3) corresponds to free- streaming (collisionally-hroadened) pre- 
equilibrium momentum space anisotropy and 5 = corresponds to thermal equilibrium. As 
in Ref. [24] , a transition width 7"^ is introduced to take into account the smooth transition 
from non-zero value of 5 to 5 = at r = Tiso- The time dependence of various quantities are, 
therefore, obtained in terms of a smeared step function [25]: 



A(r) = i(tanh[7(T - Tiso)/Tiso] + 1). 



(29) 



For r << Tiso(>> Tiso) we have A = 0(1) which corresponds to free streaming (hydrody- 
namics). With this, the time dependence of relevant quantities are as follows [24]: 



^ ^ 5(1-A(r)) 
Phard(r) = T,W"{t), 



1, 



(30) 



where. 



U{t) 
U 



U{t) ^ 
U{n) ' 



3A(r)/4 X l-5(l-A(r))/2 



' ISO 

T 



(31) 



Ti is the initial temperature of the plasma and is the velocity of sound. For isotropic 
case, we have Phard — T,Tiso — so that A = 1, U{t) = Tj/r, and U{Ti) = 1. By using 
= 1/3 we recover the Bjorken cooling law [26]. We assume here that the plasma expands 
longitudinally as the effect of transverse expansion at the very early stage might be neglected. 
This has been ilhistrated in case of photon production in Ref. [27]. Since the momentum 
space anisotropy is an early stage phenomenon, this assumption is justified. Even if the 
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PT(GeV) p^(GeV) 

Figure 2: (Color Online) Nuclear modification factor at RHIC energies. The initial conditions 
are taken as Tj=440 MeV and rj=0.147 fm/c. Left panel represents the results in the full pt 
range whereas the in the right panel the results are given for 4 < < 8 GeV 



transverse expansion is important in the very early stage, it will have two effects on the 
parton energy loss : (i) The expanding geometry will increase the duration of propagation, 
and (ii) the same expansion will cause the parton density to fall along its path. These two 
effects partially compensate each other and the energy loss is almost the same as in the case 
without the transverse expansion [28]. 

Because the colliding nuclei do have a transverse density profile, we assume that the 
initial temperature profile is given by [29] 



T,;l 



2(1- r^Rl 



1/4 



(32) 



Using Eqs.(30) and (32) we obtain the profile of the hard momentum scale as 



2 (l - rVi?^ 



(33) 



In the present work it is assumed that an isotropic QGP is formed at an initial tempera- 
ture Tj and initial time Tj. Subsequent rapid longitudinal expansion leads to an anisotropic 
QGP which lasts till tiso- Now in order to estimate the initial temperature we proceed as 
follows. In case of isentropic expansion the experimentally measured hadron multiplicity can 
be related to the initial temperature and thermalization time by the following equation [30]: 



( {b„ 



(34) 
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where (dN/dy{hjn)^ is the hadron (predominantly pions) multiphcity for a given centrahty 
class with maximum impact parameter ft^j Rt is the transverse dimension of the system, Tj 
is the initial thermalization time, ^(3) is the Riemann zeta function and = (vr^/OO) (7^ is 
the degeneracy of the system created, where gk = (7/8 • 2 • 2 • Np ■ Nc + 2 ■ 8), and Nc being 
the number of colors. The hadron multiphcity resulting from Au + Au colhsions is related 
to that from pp collision at a given impact parameter and collision energy by 



/^(U)= [(l-^)(^pa.t(U)/2 + ^(A^co«(U)J^ (35) 



dNr,. 



where x is the fraction of hard collisions. {Npart) is the average number of participants and 
{Ncoii) is the average number of collisions evaluated by using Glauber model. dN^/dy = 
2.5 — 0.25ln{s) + 0.023ln^s is the multiplicity of the produced hadrons in pp colhsions at 
centre of mass energy, ^/s [31]. We have assumed that 20% hard (i.e. x — 0.20 )and 
80% soft collisions are responsible for initial entropy production. This gives the desired 
multiplicity measured at RHIC energies. For - 10% centrality (relevant for our case) we 
obtain = 440 (350) MeV for n = 0.147 (0.24) fm/c. 

We also calculate the nuclear modification factor at LHC energies at ^/s — 2.76 TeV 
and ^/s — 5.5 TeV. For the former case, the initial conditions are estimated as follows. The 
measured charged particle multiplicity density at this energy is [32] dN^h/dy = 1600 in most 
central collisions corresponding to (A^part) ~ 382 and (Ncow) ~ 1700 [33]. Using Eq.(34) we 
obtain Tj ~ 650 MeV with Tj = 1 fm/c. For LHC energies at ^/s = 5.5, the multiplicity in 
p — p collision is not known yet. We use the empirical formula for dNpp/ dy ~ 0.8 In y/s which 
yields 71=830 for n = 0.08 fm/c. 



3 Results 

For the quantitative estimates of the fractional energy loss in an anisotropic media, we, first 
consider a plasma at a temperature T = 200 MeV with effective number of degrees of freedom 
Np=2.5 with the strong coupling as constant q;s=0.3. The fractional energy loss for non- 
zero $, (^=0.5, 1, 3, 5) for hght fiavour is shown in Fig. (1) when the quark propagates along 
the direction of anisotropy. It is observed that as the anisotropy parameter (^) increases, 
the fractional energy loss increases. The enhancement factor can be better understood by 
looking at the right panel of Fig. (1) where we have plotted the ratio of the fractional energy 
loss in anisotropic media to that in isotropic case. It is seen that at low momentum the 
enhancement is more and after that it saturates for all the values of ^ considered here. For 
larger values of the anisotropic parameter, the ratio is seen to increase reaching a maximum 
value of the order of 2 corresponding to ,^ = 5. It is to be noted that the fractional energy loss 
decreases when the direction of propagation is not aligned with the anisotropy direction [10]. 

It is observed that the energy loss increases with the anisotropy parameter ^. Mathemat- 
ically, this can be understood as follows. The energy loss is proportional to the square of the 
two-body potential v{q, ^) which has a factor R{^) in the denominator. Now this quantity de- 
creases with ^ and hence the fractional energy loss increases with ^. It is also to be noted that 
the two-body potential is stronger when the parton propagates in the anisotropy direction 
and it decreases away from the anisotropy direction leading to less energy loss [10, 16, 34]. 
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Figure 3: (Color Online) Same as Fig. (2) at = 350 MeV and = 0.24 fm/c. 

Now let us turn to the calculation of nuclear modification factor at RHIC energies. For 
the time evolution of the system Eqs.(31) and (34) are used for free-streaming interpolating 
model {S = 2). The initial conditions are taken as T, = 440 MeV and rj=0.147 fm/c. The 
results for _Raa for various values of the isotropization time ,riso have been compared with the 
PHENIX data [35] in Fig. (2). It is quite clear from right panel of Fig. (2) that the value of 
Raa for anisotropic media is lower than that for the isotropic media as the energy loss in the 
former is higher by a factor of 1.2 - 2 (see Fig. (1)). It is also observed that as Tiso increases 
the value of Raa decreases compared to its isotropic value. This is because the hard scale 
Phard dccrcascs slowly as compared to the isotropic case, i. e. the cooling is slow. Also we 
have checked that as rjso increases the rate of cooling becomes slower leading to larger value 
of the energy loss. For reasonable choices of T[so, the experimental data is well described. It 
is seen that increasing the value of t^so beyond 1.5 fm/c grossly underpredict the data. We 
find that the extracted value of isotropization time lies in the range 0.5 < Tiso < 1-5 fm/c. 
This is in agreement with the earlier finding of Tiso using PHENIX photon data [36]. 

In order to see the sensitivity on the initial conditions we now consider another set of 
initial conditions, Tj = 350 MeV and = 0.24 fm/c. The result is shown in Fig. (3). In 
this case also the data is well reproduced for the values of Tiso considered. It is observed that 
to reproduce the data larger value of Tiso is needed as compared to the case of higher initial 
temperature. We extract an upper limit of Tiso = 2 fm/c in this case. 

We also compare our results to the recent measurement of Raa by the ALICE collab- 
oration at ^/s = 2.76 TeV [33]. The result is shown in Fig. (4). Note that LHC data, at 
this energy, show quite different behavior at high (increasing trend) for which RHIC data 
do not exit. The reason might be due to larger volume and larger density. However, unlike 
the RHIC data, our present model is unable to explain the LHC data (only 2-3 points can 
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Figure 4: (Color Online) Same as Fig. (2) at y/s = 2.76 TeV corresponding to Tj = 650 MeV 
and Tj = 0.1 fm/c. 



be explained with the values of T[so considered here). Thus, within the ambit of the model 
used here, it is not possible to extract Tiso accurately at LHC energies. It is therefore, nec- 
essary to modify the present approach. The possible modifications might be the inclusions 
of coUisional energy loss, path length distribution. 

From the figures for Raa (both for RHIC and LHC energies), it is observed that the 
nuclear modification factor shows some kind of saturation behavior with Tiso- We have 
checked that this occurs until Tiso = 1-5 and 2 fm/c for two sets of initial conditions used 
at RHIC. Further increase of Tiso leads to more suppression thereby underpredicting the 
data (see Figs. (2)) . Similar observations have been noted while calculating photons from 
anisotropic media [36]. 

Although for RHIC data we could extract tiso, the extraction of tiso is uncertain by 
comparing with the LHC data at y/s = 2.76 TeV. However, for the sake of completeness, we 
predict the nuclear modification factors at higher LHC energies, i. e. at y/s = 5.5 TeV. Since 
the final multiplicity at LHC energies is higher the initial temperature is higher compared to 
RHIC energies. Also the life time of the plasma is large. Moreover, the length traversed by 
the jets in the plasma is higher. As a result, one would expect that the energy loss would be 
more. Consequently, Raa will be lower for the same values of isotropization time, Tiso- The 
results for various Tiso is delineated in Fig (5). As discussed earlier, we find that the values 
of the nuclear modification factors are less as compared to the case for RHIC. 

So far we have used the interpolating models which assumes fixed initial conditions. How- 
ever, enforcing fixed initial condition results in generation of particle number and enhance 
the final multiplicity during the transition from 5 = 2 to zero. Because most of the exper- 
imental results correspond to fixed final multiplicity, we should device a mechanism which 
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Figure 5: (Color Online) Same as Fig. (3) at LHC energies with Tj = 830 MeV and = 0.08 



enforces fixed final multiplicity. To do this the initial conditions have to be varied with the 
isotropization time, i. e. one must lower the initial "temperature' for finite riso- To ensure 
fixed final multiplicity in this model one has to redefine U{t) in Eq. (31) as in [24]: 



This redefinition corresponds to a lower initial hard momentum scale {phard{Ti) < Ti) for 
Tiso > Ti- Larger value of isotropization time corresponds to lower initial hard momentum 
scale. As we shall wee, this has important consequences on the values of Raa- 

The result for fixed final multtiplicity condition at RHIC energies has been displayed in 
Fig. (6). It is observed that the value of Raa increases with tiso compared to the isotropic 
case. The reason for this is that in FFM the larger the value of Tiso is, the lower is the 
initial hard momentum scale resulting less energy loss. This gives rise to higher value of Raa 
as comapared to the fixed initial condition (FIC) well as isotropic case. It is also 

seen that for FFM isotropic value of the nuclear modification factor is closer to the data in 
contrast to the observation made in case of FIC where it is possible to extract the value of 
the isotropization time at least for RHIC energies. Thus, use of FFM does not lead to any 
firm conclusion about the extraction of Tiso- We have checked that increasing the value of 
Tiso beyond 1.5 fm/c Raa increases grossly underpredicting the data. 

We reperform our calculation corresponding to Fig. (3) for FFM and the result is shown 
in Fig. (7). Similar behaviour (as in Fig. (6)) is observed for the reason mentioned above. 
Similarly, the calculations at LHC energies can be performed using FFM model. It is easy 
to guess that Raa will be larger as compared to the case with FIC model. 



fm/c. 
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Figure 6: (Color Online) Same as Fig. (2) for fixed final multiplicity. 

4 Summary 

In this work, we have calculated the fractional energy loss due to gluon radiation in an 
infinite size anisotropic media treating the scatterer as providing a screened coulomb-like 
potential. It is shown that the presence of initial state anisotropy increases the radiative 
energy loss by a factor of 1.2 - 2 depending upon the values of the anisotropy parameter ^ 
and the momentum of the jet when the parton propagates along the direction of anisotropy. 
We then calculate the hadron pt spectrum including this effect. A phenomenological model 
for the space time evolution of the anisotropic media has been used to obtain the time 
dependence of the anisotropy parameter and the hard momentum scale. The results for 
the nuclear modification factor for various values of the isotropization time T[so are then 
compared with the PHENIX data. It is found that for FIC the data is well reproduced if 
Tiso lies in the range 0.5 < Tiso < 1.5 fm/c which is in agreement with the previous findings. 
Increasing Tiso beyond this range seems to underpredict the data. However, for FFM due to 
the lowering of the initial hard momentum scale with Tiso, Raa increases as compared to the 
case when FIC is used. By comparing with the data it is concluded that it is very difficult 
to infer which phase (isotropic or anisotropic) is favourable. 

To cover the uncertainties in the initial conditions we also consider another set of initial 
conditions corresponding to RHIC energies. It is observed that for lower initial temperature 
the upper limit of the extracted value of Tiso is slightly higher as compared to the case where 
larger value of initial temperature is used. It is well known that the velocity of sound greatly 
influences the expansion dynamics. It has been seen that the results are extremely sensitive 
to the velocity of sound. 

After fixing the model parameters from RHIC data, we calculate the nuclear modification 
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Figure 7: (Color Online) Same as Fig. (3) for fixed final multiplicity. 

factor at LHC energies both at ^/s = 2.76 and 5.5 TeV. The data at y/s = 2.76 is compared 
with our calculation. It is to be noted that the data at y/s = 2.76 show quite different 
behavior (decreasing between = 5 — 10 GeV and again increasing) than at RHIC energies. 
It is shown that the present model is unable to predict the isotropization time at these 
energies. In order to do that further refinement of the model is necessary. Note that as the 
initial temperature is higher, the nuclear modification factors are lower as compared to the 
RHIC case for the same set of parameters. 

It is to be noted that we have not included the coUisional energy loss in the present work 
as its contribution will be subleading. However, a complete calculation must include both 
the energy loss mechanisms in order to estimate riso i- e. better constrain can be imposed 
on Tiso in such case. 

It is worthwhile to mention that the path length fluctuation is an important phenomenon 
as shown in Ref. [37] while explaining the non-photonic single electron data. In addition to 
the anisotropy parameter ^, the energy loss of parton in anisotropic media also depends on 
the direction of propagation with respect to the anisotropy axis. Therefore, consideration of 
fluctuating path length might lead to additional direction dependence. In addition, use of 
path length distribution (instead of constant value) leads to surface emission effects which 
plays significant role in determining Raa- 

It is to be mentioned that we have used the average energy loss to calculate the nuclear 
modification factor. However, to improve our calculation statistical treatment of the energy 
loss has to be incorporated. In our calculation of energy loss Landau-Pomeranchuk effect 
has been neglected which is worth investigating. 

We also do not consider the recoil of the scatterer in this work. However, this condition 
can be relaxed by incorporating the recoil corrections which plays an important role as shown 
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in Ref. [18]. Also this will lead to the unstable modes which greatly influences the transport 
coefficients, such as the radiative energy loss. This will be included in future publication. 
Furthermore, the finite size effect to the radiative energy loss in anisotropic media would 
also be interesting to study. We leave these issues for now. 
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